The listing of claims will replace all prior versions, and listings, of claims in the appUcation: 



Listing of Claims: 

1 . (currently amended) A method of obtaining magnetic resonance signals with signal 
separation for at least two chemical species in a heterogeneous magnetic field using rapid 
gradient echo imaging , comprising the steps of: 

a) obtaining first magnetic resonance signals from pixels in an object having at least two 
chemical species using a first repetition time and a first echo time, 

b) obtaining at least second and third magnetic resonant signals from the pixels using 
second and third echo times, 

c) determining a signal estimate for each species and for each pixel by combining all 
signals for the pixel using a linear least squares fitting of the signals from each pixel to 
decompose the chemical species, assuming a first value of field heterogeneity {y/o)-> 

d) calculating a first error to the field heterogeneity, 

e) repeating step c) using the first value of field heterogeneity and the error from step d), 

f) repeating step d) to calculate a second error to the field heterogeneity, and 

g) updating the value of field heterogeneity and repeating steps c) and d) until an 
acceptable error is calculated. 

2. (original) The method as defined by claim 1 wherein M chemical species are present and 
step b) includes obtaining at least M+1 magnetic resonance signals for each pixel. 

3. (original) The method as defined by claiml wherein fat and water are two chemical 
species and step b) includes obtaining three magnetic resonance signals. 

4. (original) The method as defined by claim 1 wherein in step b) the magnetic resonance 
signals at time, n=l to N, for species, j =1 to M, having real, R, and imaginary, I, parts is given 
by: 

M . M 

and a least squares fitting of all signals is given by: 

p = (A^A)-^A^S 
where A is a known matrix for M species. 
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5. (original) The method as defined by claim 4 wherein: 
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where =cos(2^/^fJ, rff = sin(2^/^fj, =27rtS- K - ^"d„ - P'fC„) and 
g[ = 271 [pl -\- p^c^- p\d^ ) are the matrix elements. 

6. (original) The method as defined by claim 4 wherein error to the field heterogeneity is 
given by: 

where B is a matrix given by: 
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where 

y = [Ml/ Ap^ Ap[ Apl Api ... A/?^ Ap'J , gl=27rt^Y.(- P^d,^- p]cj^ ^nd 

7=1 

M 

7. (original) The method as defined by claim 6 wherein step a) includes obtaining signals 
from a single coil. 
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8. (original) The method as defined by claim 6 wherein step a) includes obtaining signals 
from a plurality of coils and steps b) through g) are performed for signals from each coil, and 
further including the step of: 

h) combining field heterogeneity as determined from signals for each coil. 

9. (original) The method as defined by claim 8 wherein field heterogeneity is determined by 
weighting contributions from each coil. 

10. (original) The method as defined by claim 9 wherein the weighting contribution from 
each coil is a function of the square of the magnitude of the image contributed by that coil. 

11. (original) The method as defined by claim 10 where for each pixel, at position r, the 
combined field map is: 




p=i 



where P coils collect P independent images. 

12. (original) The method as defined by claim 1 1 wherein the combined field heterogeneity 
from step h) is smoothed by passing through a low pass filter. 

13. (original) The method as defined by claim 11 wherein M images of each species are 
obtained using signals from each of M coils and the combined field heterogeneity, and then 
combining the M images using a square root of the sum of the squares of the images. 

14. (original) The method as defined by claim 1 wherein step a) includes obtaining signals 
from a plurality of coils and steps b) through g) are performed for signals from each coil, and 
further including the step of: 

h) combining field heterogeneity as determined from signals for each coil. 

15. (original) The method as defined by claim 14 wherein field heterogeneity is determined 
by weighting contributions from each coil. 



4 



16. (original) The method as defined by claim 15 wherein the weighting contribution from 
each coil is a function of the square of the magnitude of the image contributed by that coil. 

17. (original) The method as defined by claim 16 where for each pixel, at position r, the 
combined field map is: 



where P coils collect P independent images. 

18. (original) The method as defined by claim 17 wherein the combined field heterogeneity 
from step h) is smoothed by passing through a low pass filter. 

19. (original) The method as defined by claim 17 wherein M images of each species is 
obtained using signals from each of M coils and the combined field heterogeneity, and then 
combining the M images using a square root of the sum of the squares of the images. 
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